Enhancement of the Naa;Co02 thermopower due to electronic correlations 
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Using the merger of local density approximation and dynamical mean field theory, we show how 
electronic correlations increase the thermopower of Nao.7Co02 by 200%. The newly revealed mech- 
anism is an asymmetric shift of (quasi) electrons and holes away from the Fermi level, concurrent 
with an asymmetry of the respective (group) velocities. Exploiting this effect in bandstructure and 
correlation engineering may lead to a substantial increase of the thermoelectric figure of merit. 

PACS numbers: 71.27.+a,71.1.Fd 



Climate change and the prospective oil peak necessi- 
tate the discovery of new green energy sources. One pos- 
sibility in this context is to convert hithero unused ex- 
cess heat into electrical energy using thermoelectrics [l|. 
For a real breakthrough and a widespread application of 
this technology however, thermoelectric figures of merit 
ZT > 3 are needed [J, |3|. Recently, some advance was 
brought about through phonon [4] and bandstructure en- 
gineering [5] so that ZT > 1 could be achieved. Pushing 
ZT much higher requires however a thermoelectric leap, 
which is most likely brought about by discovering a new 
class of materials - at least in the case of (high-T) super- 
conductivity scientific progress went along this line with 
the discovery of novel correlated materials. At commer- 
cially available ZT < 1 , thermoelectrics are already ap- 
plied in niche markets such as radioisotope power systems 
for satellites [6] where reliability is more important than 
efficiency; and the car industry is taking encouraging first 
steps to use thermoelectrical generators 0- The materi- 
als presently used in industry are mainly Te-based and 
have been investigated most thoroughly [8] with respect 
to optimizing their thermoelectric properties, including 
the aforementioned phonon [4] and bandstructure engi- 
neering [5]. 

The good thermoelectric properties of Te-based mate- 
rials are related to their physics of slightly doped semi- 
conductors with a low effective mass and high mobility 
of the charge carriers [I, A possible improvement are 
compounds with a pudding-mold type of bandstructure 
as those of Naa;Co02 [9.] and LiRhoOziflQl - two oxides 




with remarkably good thermoelectric properties |lll-ll4|. 

Another route to good thermoelectric properties is 
through electronic correlations [TBI, [l6| which can renor- 
malize bands most substantially, so that narrow quasi- 
particle bands with a high density of states emerge. 
The renormalization factor is most dramatic in heavy 
Fermion systems which have particularly narrow Kondo 
resonances, modeled e.g. by a periodic Anderson model 
[itI . While one might expect the narrowest resonances to 
be most suitable for thermoelectrics, there is a trade-off 
since increasing the temperature T above the Kondo tem- 
perature (essentially the resonance width) will eventually 
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FIG. 1: (Color online) Scheme of a NaxCo02 thermoelectric 
module and its crystal structure. The charge carriers are con- 
fined to the hexagonal layers of Co atoms, which are inequiv- 
alent because of the (disorderly distributed) Na ions above 
and below. 

If we apply a heat gradient, electrons and holes alike are dif- 
fusing from the hot to the cold side resulting typically in a 
very small net current. For a good thermoelectric, we need 
however a large net current, i.e. a huge asymmetry between 
electron and hole diffusion. 

Due to the pudding- mold bandstructure in NaccCo02, the 
group velocity of the holes is already considerably larger than 
that of the electrons, as indicated by the larger orange (dark 
gray) arrow. Electronic correlations additionally shift elec- 
trons away from the Fermi level, as visualized by the fewer 
blue electrons and the narrower width of the blue (light gray) 
arrow. As both effects add up, the electron is much smaller 
than the hole current leading to the net current displayed and 
hence, to the good thermoelectric properties of NaccCo02. 



smear out the resonance itself. 

In this paper, we demonstrate by means of Na2;Co02, 
which already has a favorable pudding-mold bandstruc- 
ture, that the complementary effect of electronic corre- 
lations can dramatically increase the thermopower, see 
Fig. [TJ Needed are inequivalent orbitals or sites, so that 
electronic correlations can change the "particle-hole bal- 
ance" of the system drastically, shifting more electrons 
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away from the Fermi level than holes. 

For treating electronic correlations and bandstruc- 
ture effects on an equal level, we employ the local 
density approximation+d yna mical mean field theory 
(LDA+DMFT) approach |18I-[2q| and calculate thermo- 
electric response functions as described i n 12 1|. Previous 
LDA+DMFT calculations for Na3,Co02 focused 
on the spectral function, whereas we naturally aim at the 
transport properties of this material. Starting point is a 
tight-binding fit to the LDA bandstructure [25] for the Co 
aig band within the hexagonal C0O2 layers, see Fig. [TJ 
Considering transport properties (e.g. the thermopower 
S), the other d-bands can be neglected at a first approx- 
imation. In fact, LDA predicts already that these bands 
are almost filled, and this result is also confirmed ex- 
perimentally, since Cg hole pockets are absent in angular 
resolved photoemission spectroscopy (ARPES) [26l. l27| . 
There are x<l Na+ ions per formula unit, which we con- 
sider to be randomly distributed with a Na atom or va- 
cancy in trigonal prismatic coordination to the Co sites. 
As it has been suggested in Ref. |23|, this gives rise to a 
disorder potential on the Co sites - since the positively 
charged Na+ attracts the Co-electrons. This disorder po- 
tential is of essential importance for the strong electronic 
correlations of the Co (i-electrons; and as in Ref. |23|, we 
take an disorder potential of Ae = 0.55 eV. Our resis- 
tivity calculation (inset of Fig. H]) a posteriori confirms 
that a disorder potential of this strength is present; oth- 
erwise the experimental value of the resistivity at room 
temperature, mainly due to the disorder (and not to the 
electron-electron scattering) could not be explained |28|] . 
Our model Hamiltonian hence reads 
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H = - t. 
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where cf^^Cia- are creation and annihilation operators of 
the aig electrons, with i the site index and a the spin 
index and tij denotes the hopping parameter from site 
i to site 
data 
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obtained from a tight-binding fit to LDA 
. For the Coulomb repulsion within the a^g 
orbitals we take a value U = 3.5 eV and for the Na dis- 
order potential A£ = 0.55eV as in [23]. As a DMFT im- 
purity solver we use quantum Monte Carlo simulations, 
where the Na-disorder is included at the level of the co- 
herent potential approximation CPA. For the calculation 
of the response functions we proceed as in j2l|, using a 
Taylor-fit to extract the low energy behavior of the self 
energy on the real axis. 

Fig. [2] shows our /c-resolved and /c-integrated spectral 
function. For example along path HAFK through the 
Brillouin zone, the LDA bandstructure s{k) (gray line 
inllfa)) vaguely resembles the form of a pudding mold. 
The advantage of the pudding-mold bandstructure is the 
highly asymmetric group velocity of the electrons and 
holes above and below the Fermi level, respectively. This 
velocity is given by the slope Vg{k) = de{k)/dk which. 
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FIG. 2: (Color online) (a) LDA+DMFT k-resolved spec- 
trum of Nao.7Co02 at 290 K around the chemical potential 
level (cj = 0), along with the tight-binding band extracted 
from the LDA calculation (light gray), (b) LDA+DMFT k- 
integrated spectra for T = 290 K and x=0.7 with and with- 
out interaction U and disorder potential Ae (dotted line: 
U=0, As=0; dashed line: U=0, A£=0.55 eV; solid line: U=3.5 
eV,As=0.55 eV). For fixed U=0, an increasing disorder po- 
tential As=0^0.55 eV leads to more symmetric but lower 
spectral weight in the highlighted interval ±8A:sT around the 
Fermi level. For fixed A£=0.55 eV, increasing U=0^3.5 eV 
recovers spectral weight on the hole side cj < whereas elec- 
trons Ct; > lose spectral weight in the relevant (highlighted) 
area. The inset shows the k-integrated spectra at the two in- 
equivalent lattice sites, with and without an adjacent Na^ 
(T=290 K, U=3.5 eV, and As=0.55 eV). 



because of the pudding mold, dramatically changes at 
the Fermi level at least in the FK direction (not so much 
along HA). The LDA density of states, on the other hand, 
is higher above the Fermi level (for electrons) than below 
(for holes), see (Fig. [2fb)). Hence, the asymmetry in Vg 
which increases the thermopower S is partially compen- 
sated by the asymmetry in the density of states. In Fig. 
[1] there would be, on the LDA level, more electrons than 
holes contributing. 

The thermoelectric figure of merit is defined as ZT = 
S'^/{pn)T^ with the temperature T, the electrical resis- 
tivity p and the thermal conductivity n = Kei ^ tZph 
where A^eZ, i^ph denote the electronic and phononic con- 
tributions, respectively. By use of Kubos formulas in 
DMFT [30], one finds p = I/Kq, 5' = -^f^ and 
hZei = {K2 — KI/Kq)/T with the correlation functions 
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dujA^{k,uj)^uj'', (1) 
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where V is the unit ceh volume, f{uj) is the Fermi-Dirac 
distribution function. Comparing Eq. ([T]) for n = and 
n = 1, it becomes clear that due to the factor uj in i^i, 
electron {uj > 0) and hole contributions {uj < 0) push S 
in opposite directions. A high electron-hole asymmetry is 
hence mandatory for getting a large IS*]. 

Before turning to the electronic correlations, let us first 
discuss the effect of disorder (dashed line of Fig. [2fb)), 
not yet switching on the interaction U. As one can see, 
the spectral weight around the Fermi energy is smeared 
out because of the disorder and the spectral weight is 
most dramatically reduced. As one might expect from 
Eq. ([T]), this is very disfavor able for the thermopower. 
Actually, it is reduced from the disorder-free value [9] to 
only S=25 /iV/K without interaction U. In Fig. [H we 
visualize this by a narrower arrow for both electron and 
hole contribution and a smaller net current. 

This radically changes if electronic correlations are 
taken into account. Their effect can be best understood 
by hands of the inset of Fig. [2l^b), which shows the k- 
integrated Co spectrum locally resolved for sites with 
and without an adjacent Na+ ion. The Coulomb inter- 
action drive the vacancy- adjacent Co sites more towards 
half-filling. Consequently, electronic correlations become 
stronger which reflects in the sharp quasiparticle reso- 
nance and Hubbardesque side bands. In contrast, the 
Na+-adjacent Co sites become almost completely filled 
and are only weakly correlated. Together, the shift of the 
Na+-adjacent Co spectrum below the Fermi energy and 
the strong renormalization of the vacancy- adjacent Co 
sites yields a higher spectral weight on the hole than on 
the electron side. This asymmetry now adds to the Vg{k) 
asymmetry, in stark contrast to the original LDA spec- 
tral weight. Note, that the velocities itself are, to first or- 
der, not renormalized by the electronic correlations since 
these are a property of the underlying LDA orbitals. A 
complete overview of the effect of disorder and correla- 
tion on the transport properties is given in Tab. HI 

To further illustrate the changes induced by corre- 
lation, we introduce a thermopower spectral density 
s{k,Lo) = -^^v^^{k)A^k,u;)^uj, which inte- 
grated over k and uj yields the thermopower S (Fig. [3]). 
The contour lines show that, upon (artificially) switching 
on the proper t/, thermopower spectral weight is taken 
away from the electron side. As one can see in Fig.[2fb) it 
is transferred further away from the Fermi level. In con- 
trast, the reduction of spectral weight on the hole side 
is much smaller and close to the F point s{k,uj) is even 
enhanced. This correlation-induced electron-hole imbal- 
ance strongly enhances the thermopower {S increases by 
200% w.r.t. the case of U=0). 

Finally, we present in Fig. |4] the calculated ther- 
mopower as a function of T for various doping; the 
corresponding resistivity is shown in the inset. Our 
LDA+DMFT results agree well with experiment for both 
quantities including the doping and T dependence. The 




FIG. 3: (Color online) Thermopower spectral density s(k,cj) 
for T = 290K, [7 = 3.5 eV and A£ = 0.55eV together with the 
tight-binding band extracted from the LDA data (gray) . Pos- 
itive values s > (red for cj < 0) contributed by holes domi- 
nate negative values s < (blue for uj > 0) contributed by elec- 
trons. The electronic contribution is more sensitive to correla- 
tion which is emphasized by the values s = ±10^/iVK~^eV~^ 
highlighted by an isocurve, comparing U — S.beV (solid) and 
U = (dashed). 



O [14] Kaurav et al. XNa=0.7 

[31] Wang et al. XNa=0.7(inset) 
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FIG. 4: (Color online) Thermopower 5* as a function of tem- 
perature for different doping XNa = 0.6, 0.7 showing good 
agreement to the experiment of Kaurav et al. [13]. The inset 
shows the corresponding result for the resistivity compared to 
experiments by Wang et al. ^if. 



large thermopower of NaajCo02 originates from the pud- 
ding mold band structure and is further enhanced by 
electronic correlations. The physical mechanism for the 
large resistivity on the other hand, is the disorder scat- 
tering. 

From our calculations we conclude that Naa;Co02 is 
far from optimal for thermoelectric properties. On the 
bandstructure side, there is still no strong electron-hole 
asymmetry for the group velocities in the FL and HA di- 
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parameters 


p(mQcm) S{/jY/K) 


(mW/cm 


K) 5'Vp(mW/K2 cm) 


Z (10-^ K) 


ZT 


= 0, = 


0.21 


54 


36.8 


13.8 


0.135 


0.039 


U = 3.5 eV, Ae = 


0.31 


66 


27.4 


13.7 


0.146 


0.042 


[/ = 0, As = 0.55 eV 


2.51 


25 


2.9 


0.2 


0.004 


0.001 


U = 3.5 eV, Ae = 0.55 eV 


1.37 


63 


6.5 


2.9 


0.040 


0.012 


experiment 


1.04f31] 


64fl3] 


4.0f32] 


3.9[13,^1] 


0.053[13, 31, 32] 0.0 15 [13, 31, 32] 



TABLE I: Thermoelectric and transport properties for single crystalline Nao.7Co02 at 290 K. For the computed figure of 
merit Z = / p{Kph + /^e), we take the experimental phononic thermal conductivity Kph = 66mW/cm K 32]. The electronic 
part of the thermal conductivity Ke is calculated with the correlation functions from Eq. ([Tj) and estimated experimentally by 
the Wiedemann- Franz law. Note that there are no free parameters and that only disorder (Ae) and correlation (U) together 
describe all experimental values reasonably well. 



rection due to the three dimensionality, see Fig.O Hence, 
the electron contribution (blue in Fig. [2fb)) still reduces 
the thermopower considerably. This can be improved by 
bandstructure engineering the pudding mold form, e.g., 
by expanding the c-axis length via intercalation. Sec- 
ond, disorder broadens the LDA substantially and dra- 
matically decreases the Seebeck coefficient as well as the 
conductivity ^ 1/p. There is hence a high potential for 
further increasing 5, e.g., by realizing a similar physical 
situation on the basis of inequivalent orbitals or super- 
structures without disorder. Last, but not least, we have 
shown the crucial role correlation can play, not only be- 
cause of the narrowing of the coherent excitations but 
also in strengthening the particle- hole imbalance. Given 
the high sensitivity of correlated systems to small changes 
of the external conditions and recent advances in engi- 
neering man-made heterostructures of correlated materi- 
als, there is plenty of room for improving this correlation 
effect further. 

We can thus conclude from our study of transport 
properties in Naa;Co02 that, if properly combined, band- 
structure and correlation effects (as those present in tran- 
sition metal d and rare earth f orbitals) can consider- 
ably improve the thermoelectric properties of materials. 
Jointly, engineering bandstructure and correlations hence 
opens a new route to further enhance S (and possibly ZT) 
up to the levels required for a widespread technological 
applications. 
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